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The cross section for the proton weak capture reaction ^\i{p,e '^is calculated with wave 
functions obtained from a number of modern, realistic high-precision interactions. To minimize 
the uncertainty in the axial two-body current operator, its matrix element has been adjusted to 
reproduce the measured Gamow-Teller matrix element of tritium /3 decay in model calculations 
using trinucleon wave functions from these interactions. A thorough analysis of the ambiguities that 
this procedure introduces in evaluating the two-body current contribution to the pp capture is given. 
Its inherent model dependence is in fact found to be very weak. The overlap integral — 0) 

for the pp capture is predicted to be in the range 7.05-7.06, including the axial two-body current 
contribution, for all interactions considered. 

21.30.-x, 21.45.-bv, 25.10.-fs, 95.30.Cq 

I. INTRODUCTION 

The proton weak capture on protons is the most fun¬ 
damental process in stellar nucleosynthesis: it is the hrst 
reaction in the pp chain converting hydrogen into he¬ 
lium, and the principal source for the production of en¬ 
ergy and neutrinos in main-sequence stars. The theoret¬ 
ical description of this hydrogen-burning reaction, whose 
cross section cannot be measured in terrestrial labora¬ 
tories, was first given by Bethe and Gritchheld |^, who 
showed that the associated rate was large enough to ac¬ 


count for the energy released by the Sun. Since then, a 
series of calculations has refined their original estimate 
by either conmuting the required wave functions more 
accurately or by using more realistic models for the 
nuclear transition operator [^|-^. We here contribute to 
this effort by providing an integrated study of these two 
aspects with emphasis on a reliable estimate of their as¬ 
sociated theoretical uncertainties. 

This paper is divided into seven sections and an ap¬ 
pendix. In Sec. II we set up the framework for the 
present study, by providing expressions for the pp fu- 
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sion cross section and the required matrix elements and 
by summarizing the current “best” values for the vari¬ 
ous coupling constants, Fermi function, etc. In Secs. Ill 
and IV we give a fairly detailed description of, respec¬ 
tively, the pp and deuteron wave functions, as obtained 
from modern (high-precision) interactions. The latter 
include, along with the short-range nuclear part, a com¬ 
plete treatment of electromagnetic effects up to order , 
a being the fine-structure constant, and accurately repro¬ 
duce the measured low-energy pp scattering parameters 
and deuteron properties. Sections V and VI deal with 
the calculation of the pp cross section in the approxima¬ 
tions, respectively, in which only the one-body or both 
the one- and two-body parts of the axial current operator 
are retained. In Sec. VI we also review the evidence, as 
obtained from an analysis of tritium (3 decay, for the ax¬ 
ial two-body components (explicit expressions for them 
are listed in the Appendix). Because of their model de¬ 
pendence, we adopt the phenomenological approach of 
adjusting the cutoff masses in the meson-nucleon vertices 
and V to A axial coupling constant so as to obtain agree¬ 
ment with the experimental value for the Gamow-Teller 
matrix element in tritium (3 decay. The question of how 
this procedure impacts the pp cross section is also exam¬ 
ined. Finally, in Sec. VII we summarize our conclusions, 
and provide our “best” value for the pp overlap integral 
at zero energy. 


II. CROSS SECTION 


The spin-averaged total cross section for the 
^H(p,e“''r'e)^II reaction can be written in the form 


cj{E) 


mlf{E)Y,\{d,M\A^\pp)\\ ( 2 . 1 ) 

x;rel,n.r. ^ 


Here Gy is the vector coupling constant for which the 
value Gv = (1.149 39 ± 0.000 65) x 10“^ GeV^ is used, 
as obtained from a recent analysis of ft values for su- 
perallowed O'*" —> O'*' transitions |^; mg is the electron 
mass and Urei,n.r. is the pp relative velocity. The process 
is induced by the axial-vector part of the weak interac¬ 
tion Hamiltonian, and consequently only even parity pp 
states contribute to the matrix element. 

The naive expression for the Fermi function f{E) is 
given by 


f{E) = ^ [6{E + Am-E,- E^)p,E,EldE,dE, 
ml J 


f'{E-\-Am) /m-e 


dx x\/— \ 


E + Am 


( 2 . 2 ) 


of the recoiling deuteron is neglected. A more refined 
treatment of the phase-space factor includes the effect of 
Coulomb focusing of the emitted wave function as 
well as radiative corrections to the cross section. The 
latter have not actually been calculated for the present 
reaction but have been estimated to be comparable to 
those obtained for neutron decay |l|]. As a result, the 
Fermi function is parametrized as 


f{E) = 0.144(1-h9.04F;), 


(2.3) 


with E expressed in MeV. At A = 0 the expression in 
Eq. (^) gives 0.148, which is ab out 3 % larger than the 
more accurate estimate from Eq. (2.3). 

The deuteron and even parity pp wave functions are 
written as 
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Co°, (2.4) 


(r) = 47ry2 ^ ^ (k) k) 

L even Ml 

X YLM.ijHCl (2.5) 

where yl^ (x) are the normalized eigenfunctions of the 
two-nucleon orbital angular momentum L, spin S', and 
total angular momentum J with projection M; and 

denote, respectively, the eigenstates of the spin 
S and isospin T with projections Ms and Mt- The 
deuteron u{r) and w{r) , and pp xl {x; k) radial wave func¬ 
tions are obtained from solutions of a Schrodinger equa¬ 
tion with nuclear and electromagnetic interactions, the 
latter including corrections from vacuum polarization, 
magnetic moment, two-photon exchange, and Darwin- 
Foldy terms. A discussion of the interactions and radial 
wave functions is given in Secs. Ill and IV below. Here, 
it suffices to say that xl{x\ k) behaves asymptotically as 

XL(r-,k) ^ cosSLFLikr) 3-SYo.6LGLikr), (2.6) 

where (5 l is the phase shift, and El and Gl are the reg¬ 
ular and irregular Coulomb functions. 

The nuclear axial current operator consists of one- and 
two-body components 



(2.7) 

where a = ± is an isospin index, and 


i 

(2.8) 


(2.9) 


The ratio of the axial to vector coupling constants, 
PA = GaIGv, is taken to be [|ll| 1.2654 ± 0.0042 by 
averaging values obtained, respectively, from the beta 
asymmetry in the decay of polarized neutrons (1.2626 ± 


where Am = 2mp — = 0.931 25 MeV {m^ and 

md are the proton and deuteron masses, respectively), 
E = k'^/mp is the c.m. incident energy, and the energy 
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0.0033) and ft{n) and /i(0+ ^ 0+), and gA = 

[2ft{0+ 0^)//i(n) - l]/3 = 1.2681 ± 0.0033 l|nl. The 

form of the axial two-body current operator depends on 
the dynamical model used to construct it. However, the 
need for such a term is based on an analysis of tritium /3 
decay. This evidence as well as the impact of the ambi¬ 
guities associated with the form of on the pp fusion 
cross section are discussed below in Sec. VI. 

Selection rules for a vector operator rest rict the sum 
over L in the initial capture state, Eq. (|2.5| ), to the val¬ 
ues L = 0 and 2. However, the L = 2 contribution is 
negligible at very low energies. Indeed, the initial S- and 
D-wave channel contributions to the matrix element of 
the dominant operator are proportional to, respec¬ 
tively, 


X 3 pOC 

dr u(r)xo(^; ^) — / dr w(r)Fo(fcr), 

Jo 

O nOO 

dr ir(r)x2(^; ^) — / dr w{r)F2{kr), 
Jo 


( 2 . 10 ) 

( 2 . 11 ) 


where the XL(j']k) radial wave functions have been re¬ 
placed with their asymptotic forms by setting dj;, ~ 0, 
which is appropriate for the energy range under consid¬ 
eration here (a few keV). It is then easily seen that 


J^drw{r)F2{kr) ^ // J_\ / 

J,°^dru{r)Fo{kr) " V V J 

j^dry/rw{r)I^ {2^ampr) 
J^dry/ru{r)Ii (2^ampr) ’ 

( 2 . 12 ) 


where rj = a/urei,n.r., Il are modified Bessel functions, 
and the asymptotic expressions, valid in the regime where 
rj » kr, have been used for the Fl 0- The ratio above 
is found to be roughly 0.00013 in the limit Wrei.n.r. ^ 0 
(corresponding to rj ^ oo). 

Finally, the dependence of A^ upon the momentum 
transfer q = — Pe — Pi/, where Pe and are the out¬ 
going lepton momenta, is ignored in Eq. (2.8), because 
of the very low energies involved. At if = 0 the kinetic 
energy available to the final state is only about 420 keV, 
and the finite momentum transfer correction to the ma¬ 
trix element of A*^^) for S'-wave capture can be estimated 
to be approximately [qraf' — (0.0042)^, where q ~ 420 
keV and ~ 2 fm is the root-mean-square radius of the 
deuteron—a tiny correction, indeed. 


with XL(r;fc) the radial wave function, nip the proton 
mass, and L the orbital angular momentum. The c.m. 
relative momentum k is given by = rupTiab/S. The 
boundary conditions for the wave function are 

Xz,(0;/c) = 0, 


XL{r-,k) ^-^FL{kr)Ci+GL{kr)C2, (3.2) 

with Fl and Gl the standard regular and irregular 
Coulomb functions 0- The potential V{r) can be di¬ 
vided into a long-range electromagnetic part Vem and a 
short-range nuclear part Vn ■ The coefficients Gi and G 2 
contain all the necessary information about the partial 
wave, which is usually expressed in terms of the phase 
shift S^: 


tan^*^ = G2G^ 


(3.3) 


In a practical calculation, the Schrodinger equation ( |3.1| ) 
is integrated out to some r, large in comparison with the 
range of the short-range (i.e., nuclear) force. The nu¬ 
merical wave function is then matched to the asymptotic 
form of Eq. (3.2), and the corresponding phase shift is 
defined to be the phase shift of the nuclear force with 
respect to Coulomb wave functions. Following the nota¬ 
tion of Ref. 1^, we have added the superscript C, for 
Coulomb. 

However, in reality the electromagnetic interaction in 
pp scattering is much more complicated than just the 
simple Coulomb interaction. This leads to some practi¬ 
cal problems in applying the scenario of integrating the 
Schrodinger equation, matching to Coulomb functions, 
and extracting the phase shift. This will be discussed 
below. 

The full interaction, up to second order in the fine- 
structure constant a « 1/137, is given by lTB-|Er 


Cem(pp) = Vci + Vc 2 + Evp + f 4 iM + Edf, ( 3 . 4 ) 


where 


Vci = G 


iFc{r) 


(3.5) 


= - 


2ml 


aa 

rrip 


(A + fc2) 


Fair) 


Fc{r) , Fc{r) 


(A -f fc2) 


(3.6) 


III. pp WAVE FUNCTION 


T/ ^ T ( \ 

Wp = -5--'vp(r), 

37r r 


(3.7) 


The low-energy pp scattering is described by the radial 
Schrodinger equation 


dr^ 


+ k^ 


L{L+1) 


rUpV (r) 


XL{r-,k) = 0, 


(3.1) 


Emm = — 


4mX'’ 


-Fs{r)cri-(T2 + 


Ft{r) 


S 


12 


a f, ,^Fis{r)^ ^ 


(3.8) 
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(3-9) 
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Here the Fc{r), Fs{r), Ft{r), and Fis{r) are functions 
representing the finite size of the nucleon charge distri¬ 
bution. In the limit of point nucleons, Fc{r) = Ft{r) = 
Fis{r) = 1, whereas Fs{r) = 47r(5^(r). Their explicit r de¬ 
pendence is given in Ref. . The various contributions 
are described as follows. 

The Coulomb potential Vci contains a well-known p^ ] 
energy dependence through a' = 2ka/(rripViah)■ How¬ 
ever, at the extreme low energies of interest to astrophys- 
ical calculations (a few keV), this energy dependence is 
negligible, and we can set a' = a for all practical pur¬ 
poses. 

The two-photon-exchange interaction Vc 2 behaves like 
1/r^, and so we immediately have the problem that, in 
principle, we have to integrate out to infinitely large dis¬ 
tances before we can match to Coulomb functions. 

The vacuum polarization potential hyp describes the 
augmentation of the photon propagator by an electron- 
positron pair. In the limit of point protons \Fr:(r) = 1], 
the vacuum polarization integral is given by ||^ 

/vp(r) = (l + F'j ("’ (3,10) 

with rUe = 0.511 MeV the electron mass. Including the 
finite-size effect, we get the more complicated expression 
as given by Bohannon and Heller where the expo¬ 
nential is replaced by 

— 2meTX —2meTX 

e ‘ ^ D [x)e 

- D^{x)+ ^pF>^{x)+ ^ip + 


solution of the potential W with respect to the solution 
with the potential V as the interaction. The well-known 
application of this procedure is the situation where V 
is the Coulomb potential and W is the Coulomb plus 
nuclear potential, and the phase shift Sn = is ob¬ 

tained by matching the numerical solution to Coulomb 
wave functions as in Eq. (3.2). 

If we are to include Vc 2 in IT, we run into the problem 
that we have to integrate out to inhnitely large distances 
before we can match to Coulomb functions. However, by 
including the point-nucleon limit of Vc 2 also in V, the 
asymptotic wave function can be expressed in terms of 
noninteger L' Coulomb functions, where L' satisfies 


L'{L' + 1) =L{L + l)-aa'. 


(3.12) 


The asymptotic behavior of the wave function is now 
given by 

xUr) ^ FLikr)Ci + GL{kr)C2, (3.13) 

with F^ikT) = FLi{kr) and similarly for GL{kr). The 
advantage is clear immediately: we only have to inte¬ 
grate out to distances large with respect to the nuclear 
interaction, which is only about 20 fm. But we have to 
bear in mind that now the phase shift is S^I'^^ 2 +fs+n^ 
where it should be understood implicitly that the su¬ 
perscript refers to the interactions in the point-nucleon 
limit, while the interaction denoted by the subscript in¬ 
cludes also the (short-range) finite-size effects. To make 
this clear explicitly, we have here separated off the finite- 
size effects by writing them symbolically as being due 
to some short-range potential Eps- In this notation, the 
phase s hift with respect to Coulomb functions, as defined 
in Eq. ([3.3|), now reads 


+ ^(3p + 3p^ -I- p^)D{x) 


c-P 


(3.11) 


where p = br, D{x) = [1 — {2meX/bY]~^, and b = 4.27 
fm~^ . As a matter of fact, the simple multiplication of 
Eq. (3.1C) with Fc{r) as an approximation to the in¬ 
clusion of finite-size effects (and which was adopted in 
Ref, fl^ ) already closely resembles the exact treatment 
( |3.11 ) where the finite-size effect is properly folded into 
the integral. 

The magnetic moment interaction Tum arises as a con¬ 
sequence of the nonvanishing value of the proton mag¬ 
netic moment, pp = 2.792 85/zo, while the Darwin-Eoldy 
term Vdf is a short-range potential, describing the finite 
size of the proton, where Fs{r) ^ 47r(5^(r) in the limit of 
point protons. 

Now that we have defined the full long-range electro¬ 
magnetic interaction, we can return to the question of 
how, in practice, to integrate the Schrodinger equation 
and extract the phase shift. We will restrict ourselves to 
S waves, and so the tensor and spin-orbit terms in the 
magnetic moment interaction vanish. It is convenient to 
define a phase shift S^, which is the phase shift of the 




- - AC1+C2 , eCl 

— ^Cl+C2+FS+Af — Oci-l-C2-|-FS-l-Af "F ^Cl+C2 


_ rCl-l-C2 , ^ 

~ ‘'Cl+C2+FS+Af PL, 


(3.14) 


where pL can be easily expressed in terms of the standard 
Coulomb phase shift ctl as 


PL = CTL' -<yL- {L' - L)tt/2. (3.15) 


The next step is to also include the vacuum polariza¬ 
tion. The case where we only have the Coulomb and 
vacuum polarization has been discussed in detail by Du¬ 
rand and Heller who derive expressions for the 
relevant asymptotic wave functions and vacuum polar¬ 
ization phase shift tl = '^ci+VP- Although the vac¬ 
uum polarization potential exhibits an exponential falloff, 
the small value of the electron mass means that the 
Schrodinger equation has to be integrated out to several 
hundred Fermi before the potential has dropped to suffi¬ 
ciently small values, and it is only then that the numer¬ 
ical solution can be properly matched to the asymptotic 
solution. 

The presence of Vc 2 considerably worsens the situa¬ 
tion. Since the 1/r^ behavior of Vc 2 is of longer range 
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than the exponential decay of Vyp, the Schrodinger equa¬ 
tion has to be integrated out to distances where Vyp 
is negligibly small as compared to Vc 2 - It is only then 
that we can match the numerical solution to the proper 
asymptotic solution and define the phase shift. Unfor¬ 
tunately, because of the slow falloff of the vacuum po¬ 
larization and the small magnitude of the two-photon- 
exchange contribution, we now have to integrate out to 
much larger distances. Even at a distance of 2000 fm the 
vacuum polarization has only dropped to about 1% of 
the two-photon exchange. 

The scenario of getting the pp wave function for a 
particular nuclear interaction Vn in the S wave in the 
presence of the full electromagnetic interaction Uem is 
now as follows. We integrate the Schrodinger equation 
out to a distance of 3000 fm, where the numerical so¬ 
lution is matched to the electromagnetic wave functions 
Fo{kr) and Go{kr). The latter are defined to be the so¬ 
lutions of the Schrodinger equation in the presence of the 
point-nucleon C1-I-C2-I-VP interaction. This procedure, 
therefore, determines the phase shift of the nuclear 
plus full electromagnetic interaction with respect to the 
point-nucleon Cl -|- C2 -|- VP interaction. It should be 
stressed that this phase shift is not the same as the phase 
shift of the nuclear interaction in the presence of only 
the Coulomb interaction ((5p_|_^). The relation between 
these electromagnetic wave functions and the standard 
Coulomb wave functions Fo(fcr) and Go{kr) is given by 


in the full electromagnetic interaction defines a rather 
complicated effective-range function ||^. On the other 
hand, the normalization (3.2) advocated here allows for 
an immediate substitution of the numerical wave function 
(as obtained from solving the Schrodinger equation) into 
the expression for A as defined by Salpeter [^, without 
having to worry about a phase shift which goes through 
zero at these extreme low energies and without having to 
define a complicated effective-range function. 


IV. DEUTERON WAVE FUNCTION 


The deuteron is the bound state of protons and neu¬ 
trons in the coupled ^Si+^Di two-nucleon system. For 
a given local NN potential U(r), the radial wave func¬ 
tions u{r) and w{r) for the deuteron S and D states, re¬ 
spectively, can be obtained from the coupled Schrodinger 
equation 

u{r) = m [Too(?')«(?’) + Vb 2 (r)w(r)] , 




- w{r) = m [V 2 o{r)u{r) + V 22 {r)w{r)] , 


(4.1) 


( \ = f + '^o) sin(po + '^o) W Po\ 

\Go J V - sin(po + t-q) cos(po + To) J \Go J 


where m is twice the reduced mass of proton and neutron, 
i.e., 


(3.16) 


with pq and Tq the two-photon-exchange and vacuum po¬ 
larization S'-wave phase shifts, respectively. The prime 
in the vacuum polarization is to indicate that this is 
the vacuum polarization phase shift in the presence of 
Vc 2 , which is slightly different from what is defined 
in Refs. @13 • Note that the numerical wave func¬ 


tion is now properly normalized as in Eq. (3.2), since 

JC ^ ^EM _|_ Pq _|_ 7-^, 

It should be pointed out that at extreme low energies 
(a few keV), is almost zero, po is of the order of a 
few times 10“"^ deg, whereas Tq rapidly drops from about 
— 10“^ deg at 10 keV to —10“® deg at 2 keV, and so 
exhibits a change of sign and goes through zero as a 
function of energy. Hence, it is not recommended to use 
the normalization as advocated by Kamionkowski and 
Bahcall in Ref. Q, i.e., 

^~^Co[Goikr)+cotSoFo(kr)], (3.17) 

with Gq the Gamow penetration factor. In their case S, 
there is no problem (although is very small and cot 
becomes very large), because they did not include Vc' 2 - 
Furthermore, with this normalization (3.17), the overlap 
integral A, defined below, requires knowledge of the pp 
scattering length Opp , where the presence of Vc '2 and Vyp 


_ 2mpm„ 

m =- - - 

mp + mn 


(4.2) 


All NN potentials applied in this study use consistently 
the latest, very accurate, values for nucleon masses [ p^ , 
namely, 


TOp = 938.272 31 MeV, (4.3) 

m„ = 939.565 63 MeV, (4.4) 

implying 

m = 938.918 52 MeV. (4.5) 


For the NN potential acting in particular partial waves, 
we have introduced the convenient shorthand notation 
Uoo(r) = (35i|U|3.5i), Uo 2 (r) = {^Si\V\^Di), etc., where 
(fpS'i) = Voi^(f) and (fpZli) = {v). The quantity 

7 = ifc is discussed below. 

The radial wave functions are properly normalized to 
unity, 


[u'^ {r) + w'^ (r)] =1. (4.6) 

The asymptotic behavior of the wave functions for large 
values of r is 



5 










u{r) ~ Ase 


7 (r) ^ 


1 + 


( 7 r) ( 7 r )2 


(4.7) 


where As and Ad are known as the asymptotic S- and 
D-sidXe normalizations, respectively. In addition, one 
defines the “D/S'-state ratio” 77 = Ad/As- 

Other deuteron parameters of interest are the 
quadrupole moment 


1 


pOO 

/ drr^w{r) VSu{r) — w{r) , (4.8) 

Jo -I 


20 


the root-mean-square or matter radius 


rd = 


^{lo ’ (4-9) 


and the H-state probability 


pOC 

Pd = drw^(r). 
Jo 


(4.10) 


Similar to scattering, the deuteron equation, Eq. (4.1), 
is solved numerically by integrating out to some large r 
(25 fm in our case) and match ing the numerical waves 
to their asymptotic forms, Eq. (47), producing Ag, Ad, 
and 7 from which the predicted deuteron binding energy 

is extracted. _ 

As mentioned, in the Schrodinger equation, Eq. ( 0 , 
the interaction between the two nucleons is represented 
by a local potential E(r), with r = r 2 — ri the relative 
displacement between nucleons 1 and 2. However, in gen¬ 
eral, the NN potential V is nonlocal, i.e., V = y(r,r'), 
where r is the distance between the two ingoing nucle¬ 
ons and r' the one between the outgoing nucleons. A 
local potential can then be written as E(r,r')|iocai = 
i5(r — r')y(r). For the more general case of a nonlocal 
potential, the coupled Schrodinger equation reads 


dr 2 


pOO 

i{r) =rn / dr' rr' [Voo(r,r')u(r') 

Jo 


+ Vo2il' 


Kr')], 


dr^ 


A' 


'j(r) = m / dr'rr'\y 2 o{r,r')u{r') 

Jo 

+V 22 {r,r')w{r')]. (4.11) 


This system of coupled integro-differential equations is 
then solved by a combination of finite-difference, integral- 
discretization, and matrix-inversion techniques. 

Alternatively, one may consider the two-nucleon sys¬ 
tem in momentum space, where the deuteron wave func¬ 
tion is given by 

(q) = [^o(g)%T(q) + V' 2 (g)J^ 2 T(q)] Co, (4-12) 


with the normalization 

pOO 

/ [V'o(9)+ '^ 2 ( 9 )] = 1- 


(4.13) 


The momentu m-spa ce Schrodinger equation that cor¬ 
responds to Eq. ( 4.11 ) consists of two coupled integral 
equations. 


7 + 9 Jo 


dq' 9 [Voo{q,q')'>Po(.q') 


+Vo2{q,q')'il)2{q')] , 


7 ^ 2 ( 9 ) = - 


m 


7" + Jo 


dq 9 '^ [V2o{q,q)'>Poiq') 


+V22iq,q')tp2iq')] ■ 


(4.14) 


Considering a finite set of discrete arguments for the func¬ 
tions on the left-hand side (LHS) and using the same 
set of momenta to discretize the integrals on the RHS 
produces a matrix equation that is solved easily by the 
matrix-inversion method 

The relevant Fourier transforms linking the 
configuration-space and the momentum-space ap¬ 
proaches are 

2 poo 

VLL'{q,q') = - / drr^ / dr'r 
Jo Jo 

X jL{qr)VLL'{r,r')jL'{q'r'), (4.15) 

with VLL'{r,r')\ioca.\ = 5{r - r')VLL'{r)/rr' if the poten¬ 
tial is local, and 

^^77 = , /7 /" dqq^jL{qr)ipL{q), (4.16) 

r V tt Jo 

with uo(r) = u(r), U2(r) = w(r), and jT the spherical 
Bessel functions. 

Since high reliability and precision is an important as¬ 
pect of our present investigation, we have calculated the 
deuteron wave functions fo r so me local potentials both 
ways: first, by so lving Eq. (4T) directly and, second, by 
solving Eq. (4.14) by matrix inver sion and then perform¬ 
ing the transformation, Eq. (4.16), numerically. We find 
agreement between the resulting deuteron waves to at 
least six significant digits for any r in the range 0.05- 
14 fm. This establishes the reliability of our numeri¬ 
cal methods. It also implies that in cas es wh ere we use 
the momentum-space approach and Eq. ( [4.16|) , as for the 
nonlocal potentials, our deuteron waves are of the highest 
numerical precision. 

The deuteron is a pole in the S matrix at k = 77 . 
The relativistic relationship between 7 and the deuteron 
binding energy Bd is given by 0 

\Js = rud = rup + rrin - Bd = _ ^2 _|_ ^ml - 7 ^, 

(4.17) 
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where rud denotes the deuteron rest mass. Notice that 
this equation determines the correct empirical 7 , since 
nature is relativistic. We note that in NN scattering 
we use the relativistic relationship between k and Tiab, 
which implies that the c.m. kinetic energy T is related to 
k according to 


^/s = mp + mn + T = + + k"^. (4.18) 


Thus, consistency with the scattering problem requires 
the use of Eq. (11.17 ) to determine 7 . The formal solu¬ 
tion of Eq. ([4.1^ is 


7 ^ = - {ml -ml- ml)^] /Ami, (4-19) 


and, using Bd = 2.224 575 MeV and he = 197.327 053 
MeV fm, the accurate numerical value for 7 comes out 
to be 


7 = 0.231538 0 fm 


-1 


(4.20) 


To obtain s ome pedagogical insight into 7 ^, one may 
rewrite Eq. (4.19) in factorized form 

4m^7^ = [(to„ -k mpf - ml] [ml - (m„ - mpf] 

= BdiAm-Bd){ml-6m'^), (4.21) 

where we introduce the average nucleon mass, 


= 938.918 97 MeV, 


(4.22) 


and the nucleon mass difference Sm = — mp = 

1.293 32 MeV, and used md = 2m — Bd- From this we 
get 


7 = mBd 1 - 


Am 


( (5m^\ 

V / 


(4.23) 


and rew riting twice the reduced nucleon mass [cf. 
Eq. (^)] in terms of the average mass 




(4.24) 


we finally obtain 


7 = "i-Sd 1 - 


Bd\ 1 - bm^/ 


TOt 


mBd 


Am /I — bm?/{Am?') 
4m \ m^ ) 


mBd 1 - 


Am 


(4.25) 


The approximations involved in Eq. (4.25) are good to 
1 part in 10®. Therefore, this equation rep rodu ces the 
exact value for 7 to all digits given in Eq. (4.20). One 
can now identify the term mBd as the nonrelativistic ap¬ 
proximation to 7 ^ and the factor (1 — Bd/Am) as the 


essential relativistic correction. In most calculations of 
the past, the nonrelativistic 7 was used, 7 nr = ^/rhBd = 
0.231606 6 fm“^. The difference between and the 
correct 7 , Eq. ( 1 ^) , leads to a small difference (0.03%) 
in the overlap integral (see below). Although the dif¬ 
ference is rather small, we believe one should use the 
relativistically correct value, Eq. ( 1.2C| ). 

Besides the strong interaction, there is also a nonvan¬ 
ishing electromagnetic interaction between protons and 
neutrons that can be written as llql 


VEM{np) = Vci{np) + Vuminp), 


(4.26) 


where 


Vciinp) = a/3, 


Fnp{r) 


Vuminp) = - 


Ampmn 




[Fs{r)cri-cr2 -b 


Ft{r) 


S 12 


« Fis{r) .T O , T A \ 

Mn-3—(L-S -b L-A). 


mpm 


(4.27) 


Here A = (cri — 0 - 2 )/2, and Fnp{r) is a short-range func¬ 
tion representing the finite size of the neutron charge dis¬ 
tribution (for details, see Ref. Hi)- Because the S'-wave 
expectation values for the tensor and spin-orbit opera¬ 
tors vanish, the long-range 1 /r^ parts do not contribute 
for L = 0. For L 7 ^ 0, we make the approximation that 
'^EM-I-N ~ (or Ri Sn in terms of the S ma¬ 

trix). This means that in our calculations the asymp¬ 
totic beh avio r of th e deu teron wave functions still satis¬ 
fies Eqs. (g)^ and (ini). 

The interaction (4.26) is included only in the case of 
the Argonne AV18 NN potential where it contributes 
18 keV to the deuteron binding energy, mostly from 
the magnetic moment part Vmm(’^ 7 ') of the interaction. 
One would expect that the inner part of the deuteron 
wave function is affected by the inclusion or omission of 
VEuinp) (the outer part is essentially insensitive since it 
is ruled by 7 which is identical for all potentials). For¬ 
tunately, it turns out that the quantitative effect is very 
small, as will be demonstrated below. Thus, also mod¬ 
els that do not include the electromagnetic interaction 
between protons and neutrons can be considered as suf¬ 
ficiently reliable for our study. 

Since all N N potentials are fitted to the value of 7 given 
in Eq. ( [4.2C| ), they all accurately describe the empirical 
deuteron binding energy, Bd = 2.2 24 575(9) MeV ^], via 
the relativistic relation Eq. ( 4.17 ). The other deuteron 
parameters, as well as the ^Si scattering length at and 
effective range r^, are listed in Table ||. Predictions 
are given for the five high-precision NN potentials that 
we focus on, namely, AV18 CD-Bonn |Q, Nijm- 
I iQ, Nijm-II iQ, and Reid93 Q. Notice that not all 
quantities in Table | are independent. For example, the 
deuteron effective range pd = p{—Bd,—Bd) is related to 
As, 77 , and 7 by 
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AUl + V^) 


27 

1 - 7Pd' 


(4.28) 


Gamow penetration factor 
conventionally defined as 


and the overlap integral is 


For our present investigation, essentially only As is of 
relevance (besides 7 ). However, As (and pd) cannot be 
measured directly. The empirical information given in 
the last column of Table | on As and pd are model- 
dependent extrapolations of low-energy data. Therefore, 
to trust the predictions for Hs by our NN potentials, 
it is important that these models reproduce accurately 
all measured low-energy data, which is confirmed by Ta¬ 
ble |. The only exceptions are the deuteron matter ra¬ 
dius rd and the quadrupole moment Qd-, which are both 
underpredicted by all potential models. There are, how¬ 
ever, meson-exchange current contributions and relativis¬ 
tic corrections for rd and Qd which may make up for the 
discrepancies |^,|3|. The D-state probability, Pd, that 
is listed in the bottom row of Table is not an observ¬ 
able. It is, however, an interesting theoretical quantity in 
studies of the nuclear force. The lower value for Po pre¬ 
dicted by CD-Bonn is a reflection of the nonlocal nature 
of this potential which is based upon relativistic meson 
field theory |^,Q. Meson-exchange Feynman diagrams 
are, in general, nonlocal expressions that are represented 
in momentum space in analytic form. 

Finally, in Fig. |^, we display the deuteron wave func¬ 
tions produced by the five NN potential models. Ma¬ 
jor differences are, again, related to whether the mod¬ 
els are local or nonlocal. While the central poten¬ 
tials of AV18, Nijm-II, and Reid93 are stricly local, 
the Nijm-I central force includes momentum-dependent 
terms which give rise to nonlocal structures in the equiv¬ 
alent configuration-space potential. This affects the 
deuteron S wave and is the reason why the u{r) gen¬ 
erated by CD-Bonn and Nijm-I are so similar (large solid 
and dashed curves in Fig. and differ from the other 
three potentials. The Nijm-I tensor potential is strictly 
local, similar to AV18, Nijm-II, and Reid93, which ex¬ 
plains why these four potentials generate very similar D 
waves. The CD-Bonn tensor potential is nonlocal. 


V. AXIAL ONE-BODY CURRENT 
CONTRIBUTION 


Using the wave functions as defined in Eqs. (2.4) and 
( p.5[ ) and ignoring the D-wave contribution in the initial 
scattering state, we find that the matrix element of the 
(dominant) one-body part of the axial current is given 

by 


gi()o 

((i,M|AW|pp) = J dru{r)xo{r-,k) 

PACoAiE), (5.1) 


327r 


where Ajj^ 2 ±i,o ^be spherical components of A^^), 
= T(Aa; ± iAy)/V2 and A^^o = A^, Cq is the 


A(A) = ( 7 ^/ 2 ) ^ dru{r)xo{r;k). (5.2) 


The constant 7 is defined in Eq. ( 4.20 ), and the wave 
function xo is normalized as in Eq. HI)- Because the 
solar fusion reaction actually occurs at energies of only 
a few keV, the phase shift Jq is extremely small, and so 
the exponential can conveniently be approximated 
by unity. Note that when we adopt the normalizat ion a s 
advocated by Kamionkowski and Bahcall |^, Eq. ( ^.17 ), 
we find 

/*00 

A{E) = j dru{r)xo{r-,k), (5.3) 

Jo 


where the scattering length app is defined as 
1 

^pp 


= lim Cp/ccotijo- 

dryn ^—^0 


(5.4) 


Equation ( p.3|) coincides with the definition of the overlap 
integral given by Kamionkowski and Bahcall ||^. How¬ 
ever, as stated in our discussion on the pp wave function, 
it is not at all trivial to calculate the correct scattering 
length Qpp when electromagnetic interactions other than 
the point-particle Coulomb interaction are present. 

In the following, we will present our results for the over¬ 
lap integral using realistic pp and deuteron wave func¬ 
tions. By realistic we mean that these wave functions 
were obtained by solving the scattering and bound-state 
equations using the recent high-precision NN potential 
models, the parameters of which were fitted to give an 
almost optimal description of the NN scattering data up 
to laboratory energies of 350 MeV (i.e., y^/data « I). 
The five NN models we consider consist of the AVI 8 Ar- 
gonne model 1^, the CD-Bonn model ||3^, two Nijmegen 
models, Nijm-I and Nijm-II [^, and a regularized up¬ 
date of the Reid soft-core potential |^. The AVIS po¬ 
tential was fitted including all finite- size effects in the 
full electromagnetic potential of Eq.(3.4), whereas the 
other four potentials used the point-particle approxima¬ 
tion, i.e., Fc{r) = Ft{r) = Fu{r) = 1 and Fsir > 0) = 0. 
Furthermore, the AVIS potential is the only model which 
includes the electromagnetic interaction ( 4.26| ) also in the 
deuteron. 

In Table H we show the r esult s for A^(Eiab) {E = 
Eiab/2) as calculated from Eq. (^). The integral was cut 
off at r = 50 fm, which is valid since beyond this distance 
the deuteron wave function has become extremely small, 
and so the contribution to the overlap integral becomes 
negligible. The results are shown for laboratory kinetic 
energies of 5, 4, 3, and 2 keV, which are extrapolated 
to define the result at zero energy. For each model we 
use the deuteron and pp scattering wave functions of that 
particular model. The dependence on the particular NN 
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model is found to be rather small. Taking the average 
over all hve models we find A^(0) = 6.975 ± 0.010. Leav¬ 
ing out the CD-Bonn model, which is quite different from 
the other models in that it is the only model with non¬ 
local tensor interactions, we find an even smaller model 
dependence with A^(0) = 6.970 ± 0.005. 

We again want to stress that these NN models were 
fitted including the full electromagnetic potential, and so 
the wave functions have to be calculated in the presence 
of this same electromagnetic interaction. Truncating it, 
for example by only including the standard Coulomb in¬ 
teraction, will modify the wave function and, hence, the 
overlap integral. In Table 111 we show the effect on A‘^{E) 
for different truncations of the electromagnetic part of 
the interaction. For the nuclear interaction we take the 
AV18 potential as an example. The other models show 
a similar trend. We consider four different truncations 
of the electromagnetic interaction, all for point-particle 
protons. The effect of Vc '2 is seen to be rather small: 
neglecting it increases A^(0) by only 0.0035, which is a 
0.05% effect. The proper inclusion of the vacuum polar¬ 
ization is much more important: neglecting it causes an 
almost 1% increase. 

Finally, for the AV18 potential we can also study 
the finite-size effects and the effect of Vem(»^p) in the 
deuteron calculation. Neglecting the finite-size effects 
underestimates A^(0) by only 0.08%, as shown in the 
table. Simply removing VEuinp) changes the binding 
energy to i3d(trunc) = —2.242 227 MeV, and hence the 
asymptotic behavior of the deuteron wave function. The 
consequence of this is that A^(0) increases by 0.03, al¬ 
most a 0.5% effect. However, if we first refit the binding 
energy [i.e., make a modified AV18 potential which does 
not include VEuinp), but which does have the proper 
asymptotic deuteron wave function], then the difference 
in A^(0) is only 0.001. Hence, the inclusion of VEuinp) 
under the restriction that the potential model correctly 
fits the experimental binding energy has only a small ef¬ 
fect on the overlap integral, as we alluded to earlier. 


VI. BEYOND THE AXIAL ONE-BODY 
CURRENT CONTRIBUTION 

In this section we review the procedure leading to the 
experimental determination of the Gamow-Teller (GT) 
matrix element in tritium (3 decay, and demonstrate the 
inability of calculations based on axial one-body currents 
and realistic wave functions from modern interactions to 
correctly predict this value. After a brief discussion of the 
axial two-body current operators, we address the issue 
of their model dependence by adopting the phenomeno¬ 
logical approach of constraining them to reproduce the 
experimental value of the ^H GT matrix element. We 
then calculate these two-body current contributions to 
the pp weak capture, examining in particular the ques¬ 
tion of how their associated uncertainties affect the pp 
cross section. 


A. Tritium /3 decay 


Evidence for the presence of axial two-body current 
contributions to weak transitions comes from the /3 de¬ 
cay of tritium. Its half-life can be expressed as 


(1 + Sjt)t — 


KIG\ 

fv{F)^ + fAg\{GTy^ 


( 6 . 1 ) 


where (5^=1.9% is the so-called outer radiative correc¬ 
tion, t is the half-life, and fv and /a are Fermi functions 
calculated by Towner, as reported by Simpson to 
have the values 2.8355 x 10“® and 2.8505 x 10^^ re¬ 
spectively. The experimental value for the combination 
KjGy is (6146.6±0.6) s, as obtained by Hardy et al. [|j. 
This value is actually 0.15% larger than that used by 
Simpson [^, (6137.2 ± 3.6) s, in his ^H /3-decay anal¬ 
ysis. Finally, (F) and (GT) denote the reduced matrix 
element of the Fermi and GT operators, which in the 
one-body limit are given by, respectively 

(F) = (3He||^T.,+ ||3H), (6.2) 


(GT) = (^He II II 3 h). (6.3) 

i 


Simpson reports the experimental value (1134.6 ± 
3.1) s for the combination (1 -|- 6fi)tfv- In order to 
extract a value for the tritium GT matrix element, it 
is necessary to calculate the Fermi matrix element. If 
the trinucleons were pure total T = 1/2, Mt = ±1/2 
states, then the Fermi matrix element would just be 
one. However, charge-symmetry breaking (GSB) and 
charge-independence breaking (GIB) and, more impor¬ 
tantly, electromagnetic effects in the nuclear interac¬ 
tion lead to a small correction. In the present study, 
such a correction is calculated using ^H and ^He wave 
functions obtained with the correlated-hyperspherical- 
harmonic (GHH) method from the AV18 two-nucleon 
interaction (including electromagnetic terms) and the Ur- 
bana UIX three-nucleon interaction [^. We find, ne¬ 
glecting isospin admixtures T > 3/2 (the probability of 
T = 3/2 components in ^He has been estimated to be 
about 0.0016%), 

(F)2 = 1 - e = 0.9987. (6.4) 

The present value for e is about twice that obtained 
by Saito et al. iQ in a (converged) Faddeev calcula¬ 
tion based on the older Argonne vu two-nucleon ]^ ] 
and Tucson-Melbourne (TM) three-nucleon |^^ interac¬ 
tions and phenomenological CSB and GIB terms con¬ 
strained to reproduce the observed mass difference in 
^H and ^He. However, the individual binding energies 
are underpredicted by this Hamiltonian model by about 
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3%. In contrast, the present AV18/UIX CHH wave func¬ 
tions reproduce the experimental binding energies of both 
systems within less than 10 keV (incidentally, the varia¬ 
tional CHH and “exact” Faddeev and Green’s func¬ 
tion Monte Carlo methods produce trinucleon bind¬ 
ing energies all within a few keV of each other). It is 
unclear at this point whether the difference in e values 
calculated here and in Ref. |^| is to be ascribed to bind¬ 
ing energy effects or to differences in the treatment of the 
electromagnetic, CSB, and CIB interactions (or both). 
We note that Simpson uses the value e = 0.0006 in line 
with the estimate of Saito et al. 

Using the measured half-life, and the values K/Gy = 
(6146.6 ±0.6) s, fA/fv = 1.005 29, (F)^ = 0.9987, and 
gA = 1.2654 ± 0.0042, the “experimental” GT matrix 
element is obtained: 


(GT) 


I exp 


= a/3(0.957± 0.003), 


(6.5) 

where the -v/S is from a Clebsch-Gordan coefficient. 

The experimental GT matrix element is compared with 
predictions from a number of modern Hamiltonians with 
various combinations of realistic two- and three-nucleon 
interactions in Table 0- We also give in Table ^ the 
calculated percent probabilities of the S-, S'-, P- and D- 
wave components in the ^H wave function |^J^ . A few 
comments are in order. First, the model Hamiltonians 
with the TM three-nucleon interaction are all designed 
to reproduce the experimental ^H binding energy in Fad¬ 
deev calculations by adjusting the cutoff mass in the TM 
force As already pointed out, the two-nucleon in¬ 

teractions employed in the present work are of high pre¬ 
cision, and produce fits to pp and np scattering data up 
to laboratory energies of 350 MeV with a per datum 
in the range 1.03-1.09. 

Second, in Table IV we also quote the results obtained 
using the relation 

(GT) ~ V3(Ps ± Po/i - Ps-/i). (6.6) 

where Ps, Pd, and Ps' are the probabilities of the S-, 
D- and S'-wave components in the ^H state. Use of such 
a relation implicitly assumes isospin symmetry—namely, 
that ^H and ^He form an isodoublet—and also ignores 
the contributio n of P-wave components. However, cor¬ 
rections to Eq. (^) appear to be very small, a few parts 
in a thousand. 


Third, the results listed in Table IV indicate that mod¬ 


ern interactions lead to predictions for the GT matrix 
element of tritium in the range -/S x (0.923 — 0.937), and 
therefore to an underestimate of the experimental value 
ranging, in relative terms, from 2.1% for GD-Bonn/TM 
to 3.7% for AV18/UIX. 


B. Axial two-body current model 

For the axial two-body current operator we use a 
slightly expanded version of the conventional tt- and p- 
meson exchange model first described by Ghemtob and 


Rho [Q. These are two-body currents associated with 
excitation of intermediate A resonances by tt and p ex¬ 
changes, the 7rp mechanism, and the contact nNN and 
pNN interactions. In the tables, these operators are de¬ 
noted, respectively, as Att, Ap, -np, irS, and pS. Ex¬ 
plicit expressions for them are listed in the Appendix for 
completeness. Here we only note that (i) the (nonlocal) 
momentum-dependent terms in the tt, p, and np opera¬ 
tors are retained in contrast to Ref. (ii) monopole 
form factors are included at the ttNN and pNN vertices 
with cutoff masses A.^. and Ap, respectively; and (iii) there 
is significant uncertainty in the leading Att and Ap con¬ 
tributions, since the A to A transition axial coupling 
is not known Q]. In the model adopted here, the lat¬ 
ter is related within the quark model to the nucleon qa, 
namely, pna = i6V2/5)gA- 

The present approach consists in using the simplest 
possible two-body operators that give an adequate de¬ 
scription of the longest-range mechanisms and of adjust¬ 
ing the cutoff masses within a given Hamiltonian model 
so as to reproduce the experimental ^H GT matrix el¬ 
ement. The contributions due to exchanges of heavier 
mesons, such as the Ai or renormalizations ef¬ 

fects, arising from A-isobar admixtures in the nuclear 
wave functions p^ , are neglected. However, in the next 
subsection it is argued that these approximations are not 
expected to impact in any significant way on the theo¬ 
retical predictions for the pp weak capture cross section 
once the two-body current model is constrained to fit the 
GT matrix element of tritium. 


C. Axial two-body current contributions to the pp 
capture and GT matrix element 


In Table VI we quote the contributions to the GT ma¬ 
trix element obtained with the CHH AV18/UIX trinu¬ 
cleon wave functions from the individual components of 
the axial current operators listed in the Appendix. The 
small differences between the present results and those 
reported in Ref. [p| are due to the slightly different val¬ 


ues used for A,r (A,r = 4.80 fm“^ in the present work 
versus A,r = 4.65 fm“^ in Ref. ||^) and, presumably to 
a lesser extent, to the fact the older calculations were 
based on a different Hamiltonian model, consisting of the 
Argonne Vi 4 two-nucleon and Urbana-VHI three-nucleon 
interactions, which, however, did reproduce the experi¬ 
mental binding energies of the trinucleons in a 34-channel 
Faddeev calculation The cumulative value for the 

calculated GT matrix element is -s/S x 0.964, about 0.7% 
larger than experiment. A slight adjustment in the cut¬ 
off masses A.^. and Ap or AA axial coupling (or both) is 
thus required to bring theory and experiment into perfect 
agreement. We will return to this point later, in Sec. VH. 

To test the model dependence, we have calculated the 
leading Att contribution with 42-channel Faddeev wave 
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functions obtained from the Hamiltonian models dis¬ 
cussed earlier, and the results are listed in Table fV. Both 
the one-body and Att contributions show a strong corre¬ 
lation with the -D-state probability in the trinucleon wave 
functions, which is obviously related to the deuteron In¬ 
state probability predicted by the underlying two-nucleon 
interaction, as is evident from Tables | and This cor¬ 
relation is a direct consequence of the dominant contri¬ 
butions due to T = 1 T = 0 ^Si-^Di (T = 0 ^Di) 

transitions for the one-body (Att) component. This has 
been verified explicitly by including only the above chan¬ 
nels in the Faddeev evaluation of the GT matrix element. 
As a result, the sum of the one-body and Att contribu¬ 
tions turns out to be essentially model independent, as 
indicated in Table 0 Such a conclusion is also expected 
to hold when the remaining two-body contributions are 
included. Thus, to reproduce the experimental GT ma¬ 
trix element, a single adjustment of the cutoff masses 
and Ap or gjvA in the axial two-body current operators 
should suffice for all Hamiltonian models considered. 

We now turn to the pp c apture. We only quote re¬ 
sults, presented in Table VH, corresponding to the AVIS 
and CD-Bonn interactions. The values calculated with 
these two models, which give the two extremes for the 
one-body contribution, 6.966 and 6.992, respectively, at 
zero energy, are within less than 0.2% when all two-body 
current contributions are included. Thus, the two-body 
part of the axial current leads to an increase of the AVIS 
and CD-Bonn one-body results, amounting, respectively, 
to 1.6% and 1.1%, consistently with the findings of the 
earlier study Q. 

Having demonstrated the model independence of the¬ 
oretical predictions for the GT matrix element and pp 
weak capture cross section, we now want to address the 
issue of how ambiguities in the axial two-body currents 
might affect this conclusion. To this end, it is useful to 
decompose the GT matrix element as 


(3He|^0..+ fo-)PnG) I'H) 

i<j 

(^He I ^ I 3 h), (6.10) 

i<j 

since the numbers of T, 5" = 0,0 and T, S' = 1,1 pairs 
in the trinucleons are much smaller than those with 
T, S = 0,1 and T, S = 1,0 [Q. It is now easy to 
see that, if (neglecting isospin-symmetry breaking cor¬ 
rections) I ^He) = Q I ^H), where Q = ti^xT 2 ,xT 3 ,x is the 
isospin-flip operator, then 

(^Hel ^Po"(g)0..+ (*j)A"(g) I 'H) 

i<j 

= f He I ^ P^{^J)0.,+ {^J)PSi^J) I 3h), (6.11) 

i<j 

since the matrix element is real, Q commutes with Pf, 
= 1, and Thus, the T = 0 and T = 1 con¬ 

tributions in Eq. (|6.7|) are expected to be of about the 
same size. This can be seen from Table where the 
sum of the T = 0 and 1 and T = 1 alone contributions to 
the GT matrix element from the individual components 
of the two-body operators are listed. 

It is interesting to define the two-body densities: 

po{x-, GT) = (^He I ^ (5(a: - nj)Ox,+ iij)Po (ij) \ ^H), 

i<j 

( 6 . 12 ) 


po{x;pp) = {pp I '^S{x-r,j)Oz^+{ij) \ d,0), (6.13) 

i<3 


such that 


(3He|^0..+ (*j) I^H) 

i<j 

= (3He|^0.,+ (*j)Pi"(*j) I'H) 

i<j 

+ (3He|^0,.+ (zj)Po"(G) I'H), (6.7) 

i<j 

where Oz,+ is the z component of any axial two-body 
current operator, and Pqi are projection operators over 
T = 0 and 1 two-nucleon states: 

PS{ij) + P[{ij) = l, ( 6 . 8 ) 

(6.9) 

In Eq. ( |6.7D most of the T = 0 (T = 1) contribu¬ 
tion is coming from conversion of a pn T, S' = 0,1 {nn 
T, S = 1,0) pair in ^H to a pp T, S = 1,0 {pn T, S = 0,1) 
pair in ^He, for example. 


/ dx po (a^) = O contribution. (6-14) 

Jo 

These densities are shown in Fig. where the po{x',pp) 
curves have been rescaled by a single factor R, obtained 
by matching the maximum of the GT and pp Att densi¬ 
ties. As can be seen from Fig. the GT and pp densities 
overlap in the region x < 2 fm. Of course, at larger 
cc-values the po( 2 :; GT) is significantly smaller than the 
po {x; pp ), O = ttS, Att, np, because of the increased bind¬ 
ing in the trinucleons. This scaling is to be expected, 
since it is a consequence of the “scaling” behavior more 
generally observed for the calculated T,S = 0,1 and 
T,S = 1,0 pair distribution functions in nuclei [Q; see 
Figs. H and Finally, we show in Fig. || the PA 7 r(a^) 
densities obtained with the AV18 and CD-Bonn Hamil¬ 
tonians for the GT and pp matrix elements. In this case, 
both the T — 0 and T = 1 contributions are included in 
the GT densities—namely, they have b een c alculated by 
removing the isospin projector in Eq. (HI)- Note that 
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the pp densities have been rescaled by a factor R ~ 39.0 
obtained by matching the maximum of the AV18 GT 
and pp densities. However, this rescaling also makes the 
CD-Bonn GT and pp densities very close (see Fig. ||), 
demonstrating that the R factor has only a very weak 
model dependence. 

The discussion above shows that two-body contribu¬ 
tions to the pp capture are essentially independent of the 
specific dynamical model adopted as long as the latter is 
constrained to reproduce the experimental value of the 
GT matrix element. 


VII. CONCLUSIONS 


We have calculated the axial matrix element for 
proton-proton weak capture using five modern high- 
precision nucleon-nucleon potentials. All these models 
give excellent fits to elastic NN scattering data with a 

/datum near 1 and reproduce measured deuteron prop¬ 
erties very well. We have paid particular attention to 
details of the electromagnetic interaction and the proper 
treatment of the low-energy pp scattering solutions. As 
noted before 0, the most important correction to the 
standard Coulomb interaction between protons is the 
vacuum polarization, which reduces the cross section by 
about 1%. We have shown that other fine details of the 
electromagnetic interaction increase the cross section by 
about 0.1%. This is in part compensated by the cor¬ 
rect relativistic treatment of the deuteron wave number 
7 , which gives a net 0.03% reduction in the cross sec¬ 
tion. Including just the axial one-body operator, the five 
models differ by only 0.3% in the calculated cross section. 

The biggest remaining uncertainty is in the contribu¬ 
tion of axial two-body currents, which can increase the 
cross section by about 1-1.5%. Three concerns were ex¬ 
pressed at the recent workshop on solar fusion rates |^] 
regarding the use of the known tritium /3-decay rate to 
predict the axial two-body current contribution to the 
pp fusion reaction: (1) the model dependence of the one- 
body contribution to the GT matrix element and the 
resulting uncertainty in the extracted two-body current 
contribution to that matrix element; (2) two-body cur¬ 
rents coupling T, = 1,0 pairs to T, = 1,1 pairs, 
which can contribute to the tritium GT matrix element 
but not to the pp capture; and (3) isobar and contact 
terms could give different contributions to the GT and 
pp-capture matrix elements, and thus knowledge of their 
sum in the GT may not be sufficient to predict their sum 
in the capture matrix element. 

Our detailed calculations show that these concerns do 
not influence the prediction of the pp-capture rate. In 
particular, (1) the model dependence in the one-body 
contribution to the GT matrix element comes mostly 
from that in the H-state probabilities. Because of the 
smaller D state predicted by the GD-Bonn potential (Ta¬ 
ble 0 ) , the corresponding prediction for this contribution 


is larger by about 1% (Table IV). However, the predic¬ 
tion obtained with this potential for the pp-capture rate 
via one-body currents is also larger by about 0.3% (Ta¬ 
ble i) because of the smaller D state in the deuteron 
(Table The axial two-body currents are necessar¬ 
ily weaker in the CD-Bonn model because they strongly 
couple the S and D states. In fact, the sum of one- 
and two-body current contributions is much less model 
dep ende nt than either as can be seen from Tables IV 
and IVIll 


(2) The axial two-body currents do not couple 
the T, = 1,0 pairs to the T, Tj, = 1,1 pairs in any 
significant way, as the discussion in the preceding sec¬ 
tion makes clear. (3) The two-body currents are large at 
small interparticle distances where nuclear forces domi¬ 
nate over binding energies. In this region the pair wave 
functions in different nuclei are similar in shape and differ 
only by a scale factor. This is the basis of the Bethe- 
Levinger conjecture Q, which can be used to relate 
processes such as pion and photon absorption, involv¬ 
ing nucleon pairs, in different nuclei Q. Thus the ratios 
of GT and pp-capture matrix elements of different two- 
body current terms are nearly the same as can be seen 
from Fig. Therefore, knowledge of their sum in the 
GT matrix element is sufficient to predict their sum in 
the pp-capture matrix element. 

Finally, as we have already mentioned, the GT matrix 
element is slightly overpredicted [\/3 x 0.964 versus the 
experimental value \/3 x (0.957 ± 0.003)]. Reducing the 
quark-model prediction for the iV to A axial coupling in 
the Att and Ap currents by 20% brings theory and ex¬ 
periment into agreement. The resulting CD-Bonn and 
AVIS values for the square of the pp overlap integral at 
zero energy are then found to be 7.045 and 7.059, respec¬ 
tively. Predictions for this quantity with other modern 
interactions are expected to fall in this range. Thus, the 
model dependence and theoretical uncertainty appear to 
be at the level of a few parts in a thousand, much smaller 
than the estimate given in Ref. [O. 
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APPENDIX A: THE AXIAL TWO-BODY CURRENT OPERATORS 


For completeness, we list here the momentum-space expressions for the axial two-body currents used in the present 
work. 

(1) Axial TT-exchange A-excitation current: 


16 


IttNN 


^ m^{rnA — m) 






X T 


Jja ' 


X kd 


' J- 


(2) Axial p-exchange A-excitation current: 


(Al) 


* (2) , ^ 4 g^il + Kp)2 /2(%) 

(q> p) 25 m2 (mA - m) 


{4 ( 


(Tj X 


kj) X kj - 


X T 


;)a (Ti X [(a-j X kj) X kj]} -I- i 


j- 

(A2) 


(3) Axial TT-exchange (pair) current: 
Ai%(q;7rS) = 


2m m2 m2 + kj 
(4) Axial p-exchange (pair) current: 

glil + KpY /2(%) 


9A kj f^{kj){{Ti X Tj)a(Ti xkj - iTj^ah + io-i X {pi+p'^)]}+ i ^ j. 


(A3) 


(2) Jp}-^] 

2 

'P ' "^3 


Aa.4(q;dS) = -gA- 


(U'.a {(cTj X kj) xkj-i [cr, X {(Tj X kj)] X (p, -k p')} 


tj y'lj g^3 ^2 j^2 \ 3 ^ ^ i ^ j ^ ) 

+ {Ti X Tj)a {qo-j • (cTj X kj) + i(crj x kj) x {pi + p') - [cr* x (cr^ x kj)] x kj}) + i^ j. 


(A4) 


(5) Axial Trp current: 


,( 2 ) 


(q; 7 rp) = -gA — 


^_ 

m (m2 


kf){ml + fc2) 


fpiki)Uikj){Ti X Tj)a [(1 -I- Kp)(T^ xki- i{pi 


(A5) 


Here q is the total momentum transfer, q = ki -|- k^, ki(j) is the momentum transfer to nucleon i (j), pi and p[ are 
the initial and final momenta of nucleon i, and /^(p)(fc)=pion (p-meson)-nucleon monopole vertex form factor. The 
quark model has been used to relate the irNA, pNA and axial NA couplings to, respectively, the ttNN, pNN, and 
gA couplings. The expression for ttS' represents the conventional pair current operator given in the literature. It is 
obtained with pseudoscalar pion-nucleon coupling. With pseudovector coupling the pion momentum kj in the first 
term in brackets would be replaced by the external momentum q and an additional term (p^ -I- p() would appear with 
the isospin structure (r^ x Tj)a- Furthermore, the pS operator includes only those terms which are proportional to 
(1 -I- KpY. Finally, m.^-, mp, m, and mA are, respectively, the pion, p-meson, nucleon, and A masses. 
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TABLE I. Triplet S-wave low-energy scattering parameters and deuteron properties. 



AV18 

CD-Bonn 

Nijm-I 

Nijm-II 

Reid93 

Empirical 

at (fm) 

5.419 

5.419 

5.418 

5.420 

5.422 

5.419(7)“ 

n = p{0, 0) (fm) 

1.753 

1.752 

1.751 

1.753 

1.755 

1.754(8)“ 

pd = p(.-Bd, -Bd) (fm) 

1.767 

1.764 

1.762 

1.764 

1.769 

1.765(5)*’ 

As (fm”R^) 

0.8850 

0.8845 

0.8841 

0.8845 

0.8853 

0.8845(9)*’ 

V 

0.0250 

0.0255 

0.0253 

0.0252 

0.0251 

0.0256(4)“ 

Vd (fm) 

1.967 

1.966 

1.967 

1.968 

1.969 

1.97535(85)'* 

Qd (fm^) 

0.270 

0.270 

0.272 

0.271 

0.270 

0.2859(3)“ 

Pd (%) 

5.76 

4.83 

5.66 

5.64 

5.70 

— 


“Reference [Pdj . 
^References ||25| | 
“Reference [|27|| . 
‘^Reference 
“References ||29|f 


TABLE 11. Square of the overlap integral A(£'iab) at vari¬ 
ous laboratory energies for the five NN potential models. The 
zero-energy results are obtained by extrapolating the preced¬ 
ing results. 


NN model 

R 

ef 

5 keV 

4 keV 

3 keV 

2 keV 

0 keV 

AV18 

1 

R 

7.002 

6.995 

6.987 

6.980 

6.965 

CD-Bonn 

1 


7.022 

7.014 

7.007 

6.999 

6.985 

Nijm I 


R 

7.002 

6.994 

6.987 

6.979 

6.965 

Nijm II 


11 

7.008 

7.000 

6.993 

6.986 

6.971 

Reid93 


11 

7.011 

7.003 

6.996 

6.989 

6.974 


TABLE III. Square of the overlap integral A(i?iab) at vari¬ 
ous laboratory energies for four different truncations of the 
electromagnetic interaction (all for point-particle protons). 
The nuclear interaction is the AV18 potential The re¬ 
sult for the full interaction with finite-size contributions is 
included for comparison. 


Vem(pp) 

5 keV 

4 keV 

3 keV 

2 keV 

0 keV 

Vci 

7.060 

7.051 

7.043 

7.035 

7.019 

Vci + Vc 2 

7.063 

7.055 

7.047 

7.039 

7.023 

Vci + Vvp 

6.993 

6.985 

6.978 

6.971 

6.956 

Vci + Vc 2 + Vvp 

6.996 

6.989 

6.981 

6.974 

6.960 

Full 

7.002 

6.995 

6.987 

6.980 

6.965 
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TABLE IV. One-body and two-body Att contributions to 
the Gamow-Teller matrix element of tritium /3 decay, obtained 
with various combinations of modern two- and three-nucleon 


interactions in CHH and 42-channel Faddeev calculations, the 


former for the AV18/UIX model only. The one-body results 
obtained from Eq. (3.6) are also quoted, while those under the 
heading “Total” give the sum of the one-body (first column) 


and Att contributions. 


Hamiltonian 

One-Body 

Eq. (1^) 

Att 

Total 

AV18 

0.924 

0.925 

0.0507 

0.975 

AV18/TM 

0.925 

0.925 

0.0546 

0.980 

AV18/UIX 

0.922 

0.923 

0.0560 

0.979 

CD-Bonn 

0.935 

0.935 

0.0427 

0.977 

CD-Bonn/TM 

0.937 

0.937 

0.0435 

0.980 

Nijm I 

0.926 

0.927 

0.0507 

0.977 

Nijm I/TM 

0.928 

0.927 

0.0534 

0.981 

Nijm H 

0.926 

0.927 

0.0504 

0.976 

Nijm H/TM 

0.927 

0.927 

0.0534 

0.981 

Reid93 

0.925 

0.926 

0.0514 

0.977 

Reid93/TM 

0.926 

0.926 

0.0549 

0.981 


TABLE V. The S-, S'-, P-, and D-state percent proba¬ 
bilities in wave functions. The results for the AV18/UIX 
model are from Ref. |^. 


Hamiltonian 

S 

S' 

P 

D 

AV18 

90.10 

1.33 

0.066 

8.51 

AV18/TM 

89.96 

1.09 

0.155 

8.80 

AV18/UIX 

89.51 

1.05 

0.130 

9.31 

CD-Bonn 

91.62 

1.34 

0.046 

6.99 

CD-Bonn/TM 

91.74 

1.21 

0.102 

6.95 

Nijm I 

90.29 

1.27 

0.066 

8.37 

Nijm I/TM 

90.25 

1.08 

0.148 

8.53 

Nijm H 

90.31 

1.27 

0.065 

8.35 

Nijm H/TM 

90.22 

1.07 

0.161 

8.54 

Reid93 

90.21 

1.28 

0.067 

8.44 

Reid93/TM 

90.09 

1.07 

0.162 

8.68 


TABLE VI. Contributions to the Gamow-Teller matrix el¬ 
ement of tritium j3 decay, obtained with the CHH AV18/UIX 
trinucleon wave functions. The cutoff masses A,r = 4.8 fm“^ 
and Ap = 6.8 fm“^ are used in the axial two-body operators. 
The cumulative result is 0.9636. The two-body results ob¬ 
tained by retaining only the contributions of the T = 1 pairs 
in tritium are also given (column labelled T = 1). 



Total 

r=i 

One-body 

0.9218 


Att 

0.0560 

0.0291 

Ap 

-0.0213 

-0.0111 

TTp 

0.0070 

0.0035 

TvS 

0.0044 

0.0025 

pS 

-0.0043 

-0.0021 
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TABLE VII. Square of the overlap integral A(i5iab) at various laboratory energies for the AV18 
and CD-Bonn interactions. The zero-energy results are obtained by linear extrapolation of those 
at Eish = 3 and 5 keV. The cutoff masses A^ = 4.8 fm“^ and Ap = 6.8 fm“^ are used in the axial 
two-body operators. The two-body contributions are added successively in the given order. 



AV18 

5 keV 

CD-Bonn 

AV18 

4 keV 

CD-Bonn 

AV18 

3 keV 

CD-Bonn 

AV18 

0 keV 

CD-Bonn 

One-body 

7.002 

7.022 

6.995 

7.014 

6.987 

7.007 

6.965 

6.985 

-|-7rS 

7.015 

7.024 

7.007 

7.016 

6.999 

7.009 

6.977 

6.987 

+pS 

7.005 

7.018 

6.997 

7.010 

6.990 

7.003 

6.967 

6.981 

d-ATT 

7.138 

7.126 

7.130 

7.118 

7.122 

7.111 

7.099 

7.089 

+Ap 

7.090 

7.092 

7.083 

7.084 

7.075 

7.077 

7.052 

7.055 

+TVp 

7.114 

7.097 

7.107 

7.089 

7.099 

7.082 

7.076 

7.060 


FIG. 1. Deuteron wave functions: large curves, u(r); small 
curves, w(r). The solid, dashed, dash-dotted, dotted, and 
long-dashed curves are generated from the CD-Bonn, Nijm-I, 
Nijm-II, Reid93, and AV18 potentials, respectively. 


FIG. 2. Gamow-Teller (solid lines) and pp (dashed lines) 
two-body densities. Note that all pp curves have been rescaled 
by a single factor, as explained in the text. 


FIG. 3. The T, 5 = 1,0 pair distribution functions for var¬ 
ious nuclei; see Ref. 1^. Note that the curves have been 
renormalized to the peak height of the density. 


FIG. 4. The T, S = 0,1 Ms = 0, ±1 pair distribution func¬ 
tions for given angles 6 between the spin-quantization axis 
and the relative position vector of the two nucleons and for 
various nuclei; see Ref. j^. Note that the curves have been 
renormalized to the peak height of the deuteron Ms = ±1 
9 — 0 density. 


FIG. 5. Gamow-Teller (solid lines) and pp (dashed lines) 
Att densities obtained with the AV18 and CD-Bonn Hamilto¬ 
nians. Note that the Gamow-Teller densities include both 
the T = 0 and T = 1 contributions—namely, they have 
been calc ulated by removing the isospin projector Pq {ij) in 
Eq. (|6.12| ). The pp densities have been rescaled by a single 
factor R ~ 39.0, obtained by matching the maximum of the 
AV18 Gamow-Teller and pp densities. 
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